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1 Introduction 

About 20% of the patients suffering from epilepsy 
cannot be treated well by anti-epileptic drugs. In 
these patients it can be a good option to resect the 
epileptic zone. In patients with brain tumors which 
are surgically treated for the tumor it might certainly 
be a good option since the epileptic focus can be 
adjacent to and/or border on the tumor tissue. 

The spatial analysis of the interictal spikes appearing 
in the MEG/EEG of these patients is an important 
tool for a clinician. This tool can localize the 
epileptic zone without the opening of the skull. 

Often spikes only constitute a very small part of the 
MEG data, making source localization complicated 
or even impossible. A better tool might be defined 
when taking other paroxysmal activity of the 
patients into the analysis such as delta activity. Delta 
activity and epilepsy seem related, e.g. Panet- 
Raymond and Gotman [1] demonstrated a left-right 
asymmetry in the delta band of the EEG in patients 
with epilepsy, and there is a lot of literature on the 
increased delta activity in these patients. Cerebral 
lesion such as cerebral tumors and ischemic strokes 
are also known to cause enhanced delta activity, 
though the whole background frequency spectrum 
may deviate from healthy controls. To clarify the 
relation between epilepsy and delta activity we study 
patients with cerebral tumors with and without 
epilepsy. The epileptic focus derived from spikes 
and the focus derived from the source localization of 
slow wave activity are compared. To complete the 
studies also fast brain waves are included in the 
analysis and examined for their possible 
significance. 

For studying the spontaneous brain activities large 
episodes of MEG data are recorded, which made it 
necessary to develop an automatic detection and 
source localization method for analyzing this data , 
similar to the method reported by Lewine and 
Orrison [2]. Only high power signal segments of 
band-passed filtered data are analyzed. 

In this manuscript we report on the implemented 
method and show some exemplary results on 
patients with cerebral tumors and healthy controls. 


2 Methods 

2.1 Data collection 

Controls were recruited from our own department 
and patients with brain tumors were recruited from 
the AZVU and gave their informed consent. In 
advance of the registration the subjects were 
instructed to relax and have their eyes closed in the 
system. Healthy controls sat upright and patients 
were lying down on a bed in the system. The MEG 
recording time ranged from 30 to 60 minutes. The 
MEG system in our laboratory is the whole head 
CTF system with 151 channels. For the protocol the 
digital sampling rate was set to 625 Hz and the anti¬ 
aliasing filter was set to 200 Hz. 

2.2 Automatic analysis 

The source analysis consists of five different parts: 

1. First the recorded signals are digitally band-pass 
filtered with cut-off frequencies according to 
standard EEG frequency bands (see Table 1). 
The analysis is performed for all six EEG bands 
for all the subjects. 


Frequency band 

Low F 

High F 

Delta 

i 

4 

Theta 

4 

8 

Alpha 

8 

13 

Low Beta 

13 

18 

High Beta 

18 

30 

Gamma 

30 

50 


Table 1: The cut-off frequencies of the standard 
EEG bands under investigation. 


2. The single moving dipole model is fit to the 
magnetic fields. Although the model with its 
low number of parameters may seem too simple 
to describe the magnetic fields accurately, its 
robustness, easy interpretability, and fast fitting 
capabilities make it an advantageous model to 
use. 





Figure 1: The distribution of data power of all 
samples (blue line) and the distribution of data 
power of samples yielding good fitting dipoles 
(red line) as a function of the data power. Data 
power is in units of the average sample power in 
delta band. 

3. To avoid analyzing all the recorded data, which 
is very time consuming, it is desirable to detect 
time samples at which the single dipole model 
has a high probability of explaining the 
magnetic field. The applied selection is based on 
an empirical relation between sample power and 
dipole error. Fig. 1 reveals this relation by 
plotting: i) the distribution of data power of all 
samples, and ii) the distribution of data power of 
samples yielding dipoles with residual errors 
below 10% (good fitting dipoles). Dipoles were 
fit to all samples of a 10 minute segment of 
MEG data, filtered between 1 and 4 Hz. In the 
MEG signals polymorphic delta activity was 
present and was recorded from a patient 
suffering from a glioblastoma. Apparently, good 
fitting dipoles correspond to samples with high 
powers. In the automatic detection and 
localization algorithm dipoles are fit only to 
samples with at least twice the average sample 
power. This appears a proper trade-off between 
the number of dipoles to fit and the resulting 
number of good fitting dipoles. 

4. MEG and MRI coordinates are matched using 
markers. In advance of the MEG registration 
three head coils are placed at anatomical 
landmarks on the head of the subject. These 
landmarks are the nasion, and pre-auricular 
points. The positions of these head coils are 
determined just before, during, and just after 
registration. The head coils are replaced by 
vitamin E capsules (visible on the MRI) and the 
MRI is scanned. With these two sets of three 
points, the MEG coordinate system is 
transformed to the MRI coordinate system. 

5. Only dipoles are visualized with residual errors 
below 10%, i.e., dipoles which explain more 
than 90% of the magnetic variance. Dipole 


density plots are created by means of voxels. 
The MRI volume is divided into 1cm 3 voxels, 
with each voxel color coded according to the 
number of dipoles it contains. Two orthogonal 
dipole projection plots summarize the spatial 
distribution of the dipoles. The dipoles are 
projected perpendicular on the MRI slices. 

3 Results 

3.1 Controls 

Though the dipole distributions among the controls 
differed somewhat, most of them share the about the 
same characteristics. Dipole distributions in controls 
are symmetrical with respect to the mid saggital 
plane, see Fig. 2. Dipoles describing delta (D) 
activity are located at posterior areas in the brain. 
The theta (T) dipoles and beta (B) dipoles are 
located close to the midpoint of the mid saggital 
plane. Dipoles describing alpha (A) activity are in 
the visual cortex area. Few beta dipoles and no 
gamma dipole fitted with residual errors below 10%. 



Figure 2: Typical dipole distributions in healthy 
controls. Delta (D) dipoles are located in the frontal 
areas, theta (T) and beta (B) dipoles are located 
close to the midpoint of the mid saggital plane. 
Alpha (A) dipoles are situated in the area of the 
visual cortex. 

3.2 Patients 

In patients the dipole distributions are often 
asymmetrical (with respect to the mid saggital 
plane). In this manuscript the results in the three 
patients A, B and C are presented. The results of 
patient A and B show most of the typical 
asymmetric dipole distributions for the different 
standard EEG frequency bands. 

Slow waves. The slow wave dipoles are in the 
hemisphere of the lesion, see Fig. 3 and Fig. 5. 











Dipoles describing delta activity are on the posterior 
margins of the lesions, and dipoles describing theta 
activity are more posterior to the lesion than the 
delta dipoles. Dipoles describing activity in the 
frequency range of the alpha rhythm are typically in 
the visual cortex just as in the case of the controls. 
However, examples of asymmetric distributions of 
dipoles in this frequency range are present. 

Fast waves. For the faster frequencies, see Fig. 4 
and Fig. 6. In the typical case the dipole 
distributions of the low beta and the high beta 
frequency band are similar and sometimes 
asymmetric. Gamma dipoles contralateral to the 
lesion have been found in several cases. 

Spikes. In the case of patient C spikes were present 
in the MEG, and were averaged and localized. In 
Fig. 7 the spikes are shown as dots and the delta 
dipoles are represented by a dipole distribution plot. 
Both dipole distributions have their maximums at 
the medial site of the lesion 



Figure 3: Dipole distributions of patient A suffering 
from a glioblastoma shown a) in a axial MRI slice, 
and b) in a coronal MRI slice. Dipoles are 
represented by a dipole projection plot. The delta 
dipoles (yellow) border on the tumor, theta dipoles 
(red) are more posterior to the tumor. Alpha dipoles 
are somewhat asymmetrical distributed in the visual 
cortex. 



Figure 4: Dipole distributions of patient A suffering 
from a glioblastoma shown a) in an axial MRI slice, 
and b) in a coronal MRI slice. Dipoles are 
represented by a dipole projection plot. The 
distributions of the low beta dipoles (white dots) and 
the high beta dipoles (purple) are asymmetric. Most 
of the gamma dipoles (green) are located 
contralateral to the lesion. 



Figure 5: Dipole distributions of a patient B 
suffering from an astrocytoma grade 1 shown a) in a 
axial MRI slice, and b) in a coronal MRI slice. 
Dipoles are represented by a dipole projection plot. 
Dipoles are colored as in Fig 3. 



Figure 6: Dipole distributions of a patient B 
suffering from an astrocytoma grade 1 shown a) in 
an axial MRI slice, and b) in a coronal MRI slice. 
Dipoles are represented by a dipole projection plot. 
Dipoles are colored as in Fig 4. 



Figure 7: Dipole density distribution of the delta 
dipoles and spikes (white dots, exact location) oj 
patient C with a cyst because of the resection of a 
glioma. The dipole density is highest medially to the 
resected area. Spikes are located at the same area. 

4 Discussion 

Performance. The presented method reduces the 30 
to 60 minutes of MEG data to simple graphs 
showing the generators of different brain rhythms. 
This is a huge reduction of data to review. The 









analysis is performed automatically on a PC or 
Unix-work station in about 3 times the recording 
time per frequency band. 

Spatial information. The routine gives plausible 
locations for dipoles: delta dipoles on the margins of 
cerebral lesions, alpha dipoles in the visual cortex, 
and gamma dipoles contralateral to the lesion. Of 
course, with a single dipole model one describes 
something like gravity points [3] of brain activity 
over curved brain surfaces and therefore not all of 
the activity described well by a single dipole gives a 
physiologically valid location. For example, beta 
dipoles are often located deep in the head at places 
like ventricles and white matter etc., and these 
dipoles probably describe widespread activity in this 
frequency band. 

Time information. With the loss of time information 
some important questions about the brain activities 
can not be answered with the simple plots: 

1. Are most of the visualized dipoles describing 
activity generated in a few bursts, or is the described 
activity continuous? 

2. Are activities in different frequencies related to 
each other. For example, are the delta dipoles time 
related with the theta dipoles? 

More advanced analysis methods are required to 
address these questions. 

Clinical relevance. The dipole distributions of 
healthy controls are symmetric with respect to the 
mid saggital plane. Therefore asymmetrical dipole 


distributions in the case of the patients seem to 
indicate that these dipoles describe pathological 
brain activity. If the spatial information of the 
dipoles has established its reliability, the 
neurosurgeon may use this information to apply a 
larger resection at the side of the epileptic focus in 
patients with brain tumors. Since epilepsy is often 
considered more a symptom of a brain pathology 
than a pathology itself, the additional information on 
the location of other pathological activity than 
spikes themselves might improve establishing a 
good clinical evaluation of the underlying brain 
pathology, and thereby improving the treatment of 
epilepsy. 
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